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The spontaneous emission rate of excitons strongly confined in quantum dots is proportional 
to the overlap integral of electron and hole envelope wave functions. A common and intuitive 
interpretation of this result is that the spontaneous emission rate is proportional to the probability 
that the electron and the hole are located at the same point or region in space, i.e. they must 
coincide spatially to recombine. Here we show that this interpretation is not correct even loosely 
speaking. By general mathematical considerations we compare the envelope wave function overlap, 
the exchange overlap integral, and the probability of electrons and holes coinciding and find that 
the frequency dependence of the envelope wave function overlap integral is very different from that 
expected from the common interpretation. We show that these theoretical considerations lead to 
predictions for measurements. We compare our qualitative predictions with recent measurements of 
the wave function overlap and find good agreement. 



I. INTRODUCTION 

Self-assembled quantum dots (QDs) are highly inter- 
esting for both applications and fundamental studies 
in many areas of optoelectronics because of their good 
optical properties and integrability with semiconductor 
nanotechnology. The latter is a significant advantage 
over atomic emitters, which implies, e.g., that the op- 
tical properties of QDs can be tuned by varying the QD 
size^"^. However, the integration in a solid-state environ- 
ment can lead to undesired effec ts suc h as dephasing^HH 
and non-radiative decay processe d 10 ^ 11 !. Furthermore, the 
physical understanding of QDs remains much inferior to 
that of atoms and thus improving the understanding of 
fundamental concepts of QDs is essential to realize the 
full potential of QD optoelectronics. 

A key optical parameter characterizing an emitter is 
the oscillator strength (OS), which describes the strength 
with which the emitter interacts with an electromagnetic 
field. For an atomic transition the OS has a fixed value, 
which is in stark contrast to QDs where the OS can be 
tuned by varying the size of the QD. Coulomb effects are 
predicted to become dominating for large QDs^, i.e. in 
the weak-confinement regime, which is relevant also for 
excitons weakly bound to impurities^. For small QDs 
the Coulomb interaction energy is negligible compared 
to the energy level spacing originating from the confine- 
ment, and in this strong-confinement regime the oscilla- 
tor strength is proportional to the square of the overlap 
integral of the envelope wave functions of the electron and 
the holeM The OS of the lowest-energy transition in QDs 
is typically larger than that of atoms by an order of mag- 
nitude and a pronounced frequency dependence of the OS 
was recently found to be due to the size-dependence of 
the envelope wave function overlap integral^ 3 -'. In this 
paper we show that very general features of the size de- 
pendence of the envelope wave function overlap integral 
lead to important predictions for the optical properties 
of QDs in the strong-confinement regime. These pre- 
dictions are confirmed by both recent experiment and 



numerical calculations using realistic parameters. Our 
results show that the common interpretation of the wave 
function overlap integral being loosely speaking equal to 
the probability of the electron and hole overlapping spa- 
tially is not correct. 



II. DEFINITION OF THE OVERLAP 
INTEGRALS 

We discuss QDs but note that our results apply for 
any quantum structure, i.e. quantum wires and wells. 
For QDs the emission energy can be tuned by varying 
the size and we shall use QD radius and optical angu- 
lar frequency interchangeably in the understanding that 
large radii lead to low emission frequencies. The OS f(ui) 
of In^G ai-a As QDs in the strong-confinement regime is 
given bj^Sl 

/M = ^JwfM, (1) 

where E p (x) = (28.8 — 7.3x) eV is the Kane energy and 
hu is the exciton transition energy. The envelope wave 
function overlap integral %f(w) is defined as 

Wu)= / d 3 rF e *(r, W )fX(r, W ) 

x J d 3 rX*(r»F e (r», (2) 

where F e (r, oS) (Fh(r,u>)) is the electron (hole) envelope 
wave function. We use normalized envelope wave func- 
tions, i.e. 

yd 3 r|F(r,w)| = l (3) 

for both electrons and holes, which implies that F(r, u) 
scales as r -3 / 2 , where r is the QD radius, which we shall 
use in the following. 



2 



Equation ([2]) is often misunderstood in that it is con- 
sidered to roughly describe the probability of measur- 
ing the electron and the hole in same point or region 
in space. It is the main purpose of this paper to show 
that this interpretation is not correct and that it leads to 
conclusions incompatible with experiments and numeri- 
cal simulations. 

We consider also the exchange integral Iex(w), which 
determines the energy splitting between dark and bright 
excitons and is of relevance when calculating the spin- flip 
rate of excitons^HH] and is given bjB^ 

7exM = / d 3 r\F e (r,u)\ 2 \F h (r,uj)\ 2 . (4) 

It is immediately clear that the exchange integral has 
the form of an overlap between the probability densities 
of the electron and the hole. Note, however, that due 
to Eq. ([3]) -Zex(w) scales as and has units of inverse QD 
volume, which means that it cannot be interpreted as 
a probability and more importantly it has direct conse- 
quences for the scaling behavior of Iex(w) as discussed 
below. 

Finally, we introduce the quantity describing the prob- 
ability that both electron and hole are measured in the 
same volume element f2, which for instance could be 
taken to be a unit cell. This quantity is given by the 
joint probability density, which takes a particularly sim- 
ple form for independent quantities^, 

ifc(w)= / d 3 v\F e (r,u)\ 2 f d 3 r'\F h (r',co)\ 2 . (5) 
Jn Jn 

We can then define the overlap probability as the prob- 
ability of measuring the electron and hole in the same 
region of space. 

P( W ) = ^Pn i (a;), (6) 

i 

where it is important to realize that -P(w) depends on 
Qi. Since the squared wave functions are probability dis- 
tributions the point probabilities vanish and thus P(u>) 
vanishes as fli becomes very small. Thus we must con- 
sider finite fii. 

III. SIZE DEPENDENCE OF THE OVERLAP 
INTEGRALS 

The size-dependence of the OS integral is commonly 
understood as follows: When the size of the QD decreases 
the electron and hole are gradually expelled from the QD 
and this is more pronounced for the electrons because of 
the lighter effective mass as compared to the heavy holes, 
which are the relevant holes because vertical confinement 
and strain lift the degeneracy with the light-hole band. 
The probability of the electron and hole coexisting at 
the same region in space is reduced and so the oscillator 
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FIG. 1. Schematic illustration of the size frequency depen- 
dence of the overlap integrals for infinite potential barriers in 
terms of the electron (E) and heavy hole (HH) wave functions. 
(a)-(c) For infinite barriers the wave functions are indepen- 
dent of the effective mass and hence follow the QD size, which 
is indicated by the black circles. The wave function amplitude 
is indicated by the color intensity. This leads to a constant 
wave function overlap (d) and characteristic size dependencies 
of the exchange overlap (e) and the overlap probability (f). 



strength must decrease for decreasing QD size, i.e. for 
increasing emission energy. Although this conclusion is 
correct the above reasoning is in fact incorrect. In the 
following we discuss this point in further detail. 

We consider the limits of very small and very large QD 
radius as well as the intermediate size regime in order 
to reconstruct the size dependence and hence frequency 
dependence of the overlap integrals Iwf( w ), Ip an d 
-P(w). For simplicity, and without loss of generality, we 
consider a spherical potential. We note that the following 
considerations are purely mathematical because, e.g., for 
large radii the strong-confinement model breaks down so 
that 7\vf(w) no longer describes the oscillator strength. 
Another issue is that for very small radii a spherical 
three-dimensional potential has no bound states while 
bound states always exist for one- dimensional potentials. 
Here we are not concerned with these issues because we 
are interested in the mathematical limits of the overlap 
integrals in order to reconstruct their frequency depen- 
dences and, as we will show, the mathematical consider- 
ations lead directly to physical predictions confirmed by 
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FIG. 2. The same illustration as Fig. [T] but for finite barri- 
ers, (a) The wave functions expelled strongly from the QD 
for very large and very small QDs. (b) In the intermediate 
regime the lighter electron mass leads to a significant frac- 
tion of the electron residing outside the QD while the hole 
is more confined, (c) For very large QDs the wave functions 
are approximately identical. This leads to the wave function 
overlap exhibiting a minimum (d) while the exchange integral 
(e) and overlap probability (f) exhibit a maximum. 



recent experiments. 

In Fig. [I] we consider infinite potential barriers, which 
are of little physical relevance but important for the un- 
derstanding. In this case the electron and hole envelope 
wave functions are identical and independent of the QD 
size as indicated in Figs, [ija)-(c). Therefore the wave 
function overlap integral is unity and independent of the 
QD size as shown in Fig. [ljd) . The exchange overlap in- 
tegral scales as inverse volume and therefore it vanishes 
for very large wave functions and diverges cubically for 
very small QD radii as shown in Fig. |TJe) . For a given 
fli the probability of measuring the electron and hole in 
the same region becomes P(uj) = 1 when the QD radius 
becomes smaller than fii while P(uj) = when the QD 
radius goes to infinity as shown in Fig. [ljf). 

Let us now examine the more realistic situation of fi- 
nite barriers. We can gain physical insight into the prob- 
lem by considering the length scales entering the system, 
namely the QD radius and the penetration length L of 
the wave functions into the surrounding material, which 



Quantum dot radius (nm) 

7 6 5 4 




1.15 1.20 1.25 1.30 

Emission energy (eV) 

FIG. 3. Measured wave function overlap (circles) for different 
emission energies compared to numerically calculated wave 
function overlap values (black curve, left axis). Using the 
same numerical parameters we have calculated the exchange 
overlap (red, right axis) and the overlap probability (blue, 
right axis), which have been normalized to their maximum 
values. 



for a one-dimensional square well is given bjE^ 

h 



L = 



(7) 



where m* is the effective mass, E is the energy, Vo is the 
confinement potential, and we have assumed Eq < Vo, 
i.e. we are considering a bound state. Here E must be 
found numerically by the wave function continuity crite- 
rion, hence L depends on energy and therefore also the 
QD radius r. We can therefore quantify three regimes. In 
the following E and HH denotes electron and heavy hole 
respectively. In the small dot regime, i.e. r <C Lhh,Le, 
the envelope wave functions are strongly expelled from 
the QD as shown in Fig.[2ja). In the intermediate regime, 
i.e. Lhh < r < Le, the electron is expelled more from 
the QD than the hole as shown in Fig. [2jb) . In the large 
QD regime, i.e. r Lhh, Le, the small fractions of the 
envelope wave functions leaking out of the QD become 
negligible and the envelope wave functions become effec- 
tively identical as shown in Fig. [5Jc). This size depen- 
dence leads directly to the size dependence of the three 
overlap integrals shown in Figs. [2]jd)-(f). In the small- 
dot regime the wave function overlap will have a finite 
value, which will be very small due to the strong depen- 
dence on the difference in effective mass for electrons and 
holes. The clear qualitative differences between the wave 
function overlap and the overlap probability in particular 
for large QD radii directly shows the incorrectness of the 
common wave function overlap interpretation. 

In order to test the findings above we have performed 
numerical calculations of the three overlap integrals for 
a QD comprised of InGaAs with 46% indium embedded 
in GaAs using a finite-element model to solve the effec- 
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tive mass equation. The result is shown in Fig. [3] The 
indium mole fraction and size distribution have been op- 
timized to fit recent e xperim ental data on the frequency 
dependence of Iwf(^V^^ as shown in Fig. [3] We note, 
however, that a similar frequency dependence of Iwf(uj) 
was found for all parameters investigated in RefP^, i.e. 
also when including various aspect ratios and strain mod- 
els. Further exp erimental and numerical details can be 
found in RefP^l. Notably, as predicted above on general 
grounds Iex(w) and P(oj) exhibit very different frequency 
dependencies from Iyyp(u)) in good agreement with the 
qualitative predictions of the previous sections. In the 
calculation of P(w) we have chosen f2j w 1 nm 3 but the 
appearance of the maximum remains for all values of fli 
that we have investigated. 

IV. CONCLUSION 

In conclusion we have shown that general mathemat- 
ical aspects of the envelope wave function overlap, the 



exchange overlap, and the overlap probability lead to dis- 
tinctly different frequency dependencies. We have con- 
firmed these predictions by numerical calculations and 
in the case of the envelope wave function overlap also 
by recent experiments. The essential conclusion is that 
the envelope wave function overlap must increase with 
increasing QD size (decreasing emission energy) which 
is even qualitatively very different from the dependence 
exhibited by the overlap probability. This shows that 
the common interpretation of the envelope wave func- 
tion overlap as an overlap probability is a misconception 
leading to predictions incompatible with experiments and 
numerical calculations. 
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